arXiv:1507.06014v2 [astro-ph.HE] 16 Dec 2015 


Draft version December 17, 2015 

Preprint typeset using DTeX style emulateapj v. 5/2/11 


SPECTRAL AND TEMPORAL PROPERTIES OE THE ULTRA-LUMINOUS X-RAY PULSAR IN M82 PROM 15 
YEARS OE CHANDRA OBSERVATIONS AND ANALYSIS OE THE PULSED EMISSION USING NUSTAR 

Murray Brightman\ Fiona Harrison^ Dominic J. Walton^’\ Felix Fuerst\ Ann Hornschemeier^’^, Andreas 
Zezas^’®’^, Matteo Bachetti®, Brian Grefenstette\ Andrew Ptak^’^, Shriharsh Tendulkar^ Mihoko Yukita^’^ 
^Cahill Center for Astrophysics, California Institute of Technology, 1216 East California Boulevard, Pasadena, CA 91125, USA 
^Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA 91109, USA 
^NASA Goddard Space Flight Center, Code 662, Greenbelt, MD 20771, USA 
"^The Johns Hopkins University, Homewood Campus, Baltimore, MD 21218, USA 
^Physics Department, University of Crete, Heraklion, Greece 
^Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 
^Foundation for Research and Technology-Hellas, 71110 Heraklion, Crete, Greece 
^INAF/Osservatorio Astronomico di Cagliari, via della Scienza 5, 1-09047 Selargius (CA), Italy 

Draft version December 17, 2015 

ABSTRACT 


The recent discovery by Bachetti et al. (2014) of a pulsar in M82 that can reach luminosities of 


up to 10^^ ergs“^, a factor of ^ 100 the Eddington luminosity for a 1.4 Mq compact object, poses 
a challenge for accretion physics. In order to better understand the nature of this source and its 
duty cycle, and in the light of several physical models that have been subsequently published, we 
conduct a spectral and temporal analysis of the 0.5-8 keV X-ray emission from this source from 15 
years of Chandra observations. We analyze 19 ACIS observations where the PSE of the pulsar is not 
contaminated by nearby sources. We fit the Chandra spectra of the pulsar with a power-law model 
and a disk black body model, subjected to interstellar absorption in M82. We carefully assess for the 
effect of pile-up in our observations, where 4 observations have a pile-up fraction > 10%, which we 
account for during spectral modeling with a convolution model. When fitted with a power-law model, 
the average photon index when the source is at high luminosity {Lx> 10^^ ergs“^) is F = 1.33 ±0.15. 
For the disk black body model, the average temperature is Tin = 3.24 ± 0.65 keV, the spectral shape 
being consistent with other luminous X-ray pulsars. We also investigated the inclusion of a soft excess 
component and spectral break, finding that the spectra are also consistent with these features common 
to luminous X-ray pulsars. In addition, we present spectral analysis from NuSTAR over the 3—50 
keV range where we have isolated the pulsed component. We find that the pulsed emission in this 
band is best fit by a power-law with a high-energy cut-off, where F = 0.6 ± 0.3 and Eq = 14^3 keV. 
While the pulsar has previously been identified as a transient, we find from our longer-baseline study 
that it has been remarkably active over the 15-year period, where for 9/19 (47%) observations that 
we analyzed, the pulsar appears to be emitting at a luminosity in excess of 10^^ ergs“^, greater than 
10 times its Eddington limit. 

Subject headings: stars: neutron - galaxies: individual (M82) - X-rays: binaries 


1. INTRODUCTION 

The discovery of coherent pulsations with a perio d of 
1.37 s in the X-ra y emission of M82 by NuSTAR (|Ba- 
chetti et al. 2014 hereafter B14), shown to be associ- 
ated with an ultra-luminous X-ray source (ULX) that 
is known to reach luminosities of 10^^ ergs“^ (with the 
assumption that the source radiates isotropically) is a 


lation as to the nature of this source (Lyutikov |2014i 

Christodoulou et al.||2014l |Eragos et al.||2U15| 

Eksi et al. 

20151 IKluzniak & Lasota|20151|Shao & Li|2015 

DairOsso 

et al. 120151 Mushtukov et al.|2015p. Since the 

pulsations 


netized neutron star with a mass of ^1.4 Mq, the ob¬ 
served peak X-ray luminosity is 100 times the system’s 
Eddington limit. 

Erom B14, the neutron star orbits its companion with 
a 2.5-day period that is close to circular, having a pro¬ 
jected semi-major axis of 22.225 light s (6.66x10^ km) 
and a companion star with a minimum mass of 5.2 Mq. 


A linear spin up is also observed from the pulsations dur¬ 
ing the NuSTAR observations, with a pulse derivative 
P —2 X 10“^^ s s“^, that varies from observation to 
observation. The pulse profile is also close to sinusoidal 
(B14). Any theoretical model must be able to account 
for all of these properties. 

High T-field {B ± 10^^ G) accreting pulsars have been 
observed to emit in excess of their isotropic Eddington lu¬ 
minosities. Magnetic fields allow pulsars to exceed their 
Eddington luminosities by tunneling the accreting mate¬ 
rial along the magnetic field lines onto the magnetic poles 
of the neutron star, with the X-ray emission radiating out 
the sides of the accretion column. Observational evidence 
that showed the pulsar SMC X-1 to be super-Eddington 
motivated calculations of radi ative transfer in the pres - 
ence of these magnetic fields by Basko & Sunyaev (|1975). 
The authors calculate the limiting luminosity of these 
systems, showing that this dep ends strongly on the ge ¬ 
ometry of the accretio n channel ( Basko fc Sunyaev|1976 ). 


DairOsso et al. (2015) consider the observational proper¬ 
ties of the pulsar in 1V182 and numerically solve the torque 
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equation with respect to the scenario in which matter is 
tunneled along the magnetic field lines, and argue in favor 
of a high magnetic field strength (B^ 10^^ G) to explain 
the variation in P and the pulsar’s high luminosity. In 
this case the strength of the magnetic field would disrupt 
the accretion disk at much larger radii (< (80 —90)i?Ns): 
causing the disk temperature to decrease. This would 
be difficult to observe, however, since the X-ray emis¬ 
sion from the accreting pulsar is dominated by the base 
of the accretion column near the neutron star’s surface. 
A magnetic field of B^ 10^ ^ G or stronger w as also fa¬ 
vored for the ULX pulsar by Eksi et al. (2015), also based 
on calculations of torque equilibrium. Magnetic fields of 
these strengths have the power to reduce the Gompton 
scattering cross se ction, thus increasing the critical ac¬ 
cretion luminosity. Dall’Osso et al. (2015) calculate that 

10^^ G field would decrease 
50. 


for 30 keV photons, a B 
the scattering cross sec tion by a factor of > 


Conversely however, Kluzniak & Lasota (2015) argue 
for a low magnetic field strength (B< 10^ G). They base 
their conclusions on the ratio of spin-up rate to lumi¬ 
nosity (10“^^ (erg s)“^), which is an order of magni¬ 
tude lower than typical X-ray pulsars. These authors 
argue that a disk truncated at large radii would not 
provide the required lever arm to power the observed 
spin-up. Another interpreta tion of the magnetic field by 


Christodoulou et al. (2014) implies that the magnetic 
field is in fact typical of pulsars at 10^^ G and that 
the observed luminosity can be accounted for by geomet¬ 
ric beaming. However, the near-sinusoidal shape of the 
pulse profile suggests that strong beaming does not oc¬ 
cur, and further, Eddington ratios of ^100 are difficult 
to reconcile with this scenario. It is clear from these ar¬ 
guments that further observational evidence, such as the 
duty cycle of the source and X-ray spectral properties 
are needed to gain more insights into the nature of the 
source and its evolution. 

The ULX associated with the pulsa r was first resolved 


by C handra HRG in October 1999 (Matsumoto et al. 


2001) and designated CXOM82 J095551.1+694045. This 


source is now kno wn to be the second most luminous X- 
ray source in M82 (lEeng & Kaaret 120071 iKong et al. 12007) 
after the ULX M8^ X-1 (GXOU ^0955^0.2+6940 47) and 
thus is c ommo nly r eferred to as M82 X-2. Both Eeng fc| 


Kaaret (2007|) and Kong et al. (2007) use the Chandra 
data on 1V182 to study this source, noting its high X-ray 
luminosity and its month-timescale variability, identify¬ 
ing it as a transient . As fo r associations at other wave¬ 
lengths, Kong et al. (2007) found that X-2 is coincident 
with the position of a star cluster seen in a near-infrared 
HST NICMOS E160W image, also associated with th e 
radio source 42.21+59.2 from McDonald et al.j ([2002[), 
i dentif ied as an Hll region. Eurthermore, [Gandhi et al^ 
(2011) presented high-resolution mid-infrared imaging of 
the center of M82 and tentatively assigned their source 
7^11 as a counterpart to X-2, based on its proximity to 
the position of the radio source. M82 X-1, which gener¬ 
ally dominates the X-ray emission of the galaxy, is a can¬ 
didate for an intermediate mass black hole, based on its 
extr eme X-ray luminosi ty, which reach up to 10^^ erg s“^ 
(e.g. Kaaret et al.pOOO ) and the detection of twin-peaked 
uasi-periodic oscillat ions at frequencies of 3.3 and 5.1 Hz 
Pasham et al. 2014). M82 X-1 and X-2 are separated 


by only 5" and thus only resolva ble by Chandra. 


Si nce the results pr esented in Eeng & Kaaret (2007) 
and Kong et al. (2007), M82 has been observed by Chan- 
dra on 18 further occasions, including the Chandra data 
used in B1 4 to identify the sourc e of the ultra-luminous 
pulsations. Eeng & Kaaret (2007) and Kong et al.j (2007) 
identify X-2 as a transient, however, B14 found that the 
source retains its high luminosity 7 years later. If X- 
2 persists at luminosities of ^ 10^^ ergs“^, assuming a 
mass to energy conversion efficiency of unity, the neutron 
star will grow at a rate of ^ 2 x 10“^ Mq yr“^, mean¬ 
ing it will collapse into a black hole within ^10 million 
years. 

The goal of this paper is to better understand the duty 
cycle of the source and its spectral characteristics, and 
discuss these characteristics with respect to theoretical 
models and other luminous pulsars. While only Chandra 
data can be used to spatially resolve the pulsar, NuS- 
TAR^ with its timing capabilities, allow it to temporally 
isolate the pulsed component due to the coherent pulsa¬ 
tions emitted by the source. We carry out spectral anal¬ 
ysis of the source in the 0.5—8 keV band using archival 
Chandra data and spectroscopy in the 3—50 keV band us¬ 
ing NuSTAR. In Section we describe the Chandra data 
and analysis and in Section we describe the NuSTAR 
data and analysis. In Section [^we present our results 
and in Section we discuss ourTindings with respect to 
other luminous X-ray pulsars and their implications for 
theoretical models. A distance of 3.3 Mpc to M82 is 
assumed throughout ( Eoley et al.||2014~ ). 


2. CHANDRA DATA AND ANALYSIS 


In total, M82 has been observed 28 times by Chan¬ 
dra^ first on 1999 September 20 with ACIS-I and most 
recently on 2015 January 20 with AGIS-S, covering more 
than 15 years of activity in the galaxy. The data con¬ 
tain a range of exposure times for each obsID, ranging 
from 2—120 ks, taken with all three instruments, ACIS- 
I, ACIS-S and HRG, all without the use of the gratings. 
Eurthermore, the galaxy has been placed at a mixture 
of on-axis and off-axis angles, the off-axis angles being 
used to mitigate pileup from the bright X-ray sources 
by spreading out the counts over a wider area of the 
detector. This has often been combined with sub-array 
readout, typically l/8th of one AGIS-S chip, also to mit¬ 
igate pileup by reducing frame times from 3.2 s to 0.4 s. 
These observational data are summarized in Table [H 

We also present images of the central 20"x20" region 
of M82 from each of the 28 obsIDs, centered on X-2 in 
Fig[l] This figure also illustrates the range of Chandra 
data available on this galaxy. It is clear from this figure 
that not all the data can be used to study X-2, as in 
many cases the PSE of X-2 is blended with the PSE of 
two nearby sources to the south (i.e. obsIDs 378, 379, 
380, 6097, 8190, and 10027). 

We carry out analysis on level 2 event data that has un¬ 
dergone standard data processing by the Chandra X-ray 
Genter (v8.4.4 for obsID 10025; v8.4.5 for 2933, 10026, 
10542, 10543, 10544, 10545, 10925, 11104, 11800, 13796; 
V8.5.1.1 for 5644, 6361, 15616; vlO.2.1 for 16580; vlO.3.1 
for 17578 and vlO.3.3 for 16023). We reprocess obsIDs 
361 and 1302 using the CIAO (v4.7, GALDB v4.6.5) script 
CHANDRA_REPRO. We extract the spectrum of the source 
for all obsIDs where the PSE of X-2 is not blended with 
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Figure 1. 20"x20'' 0.5—8 keV Chandra images of the central region of M82, consisting of 28 obsIDs taken over the 15 year-period 1999 to 
2015. The obslDs are written in the top left of each image. Images are centered on the position of the ultra-luminous pulsar, X-2, which is 
marked with white cross-hairs. The data were taken with a variety of instruments, exposure times, off-axis angles and detector sub-arrays, 
details of which are given in Table From top-left to bottom-right, the images are ordered by epoch, starting in Sep tember 1999 and 
ending in January 2015. The green circle to the north-east of X-2 marks the radio kinematic center of the galaxy from [Weliachew et al.| 
( |1984| . North is to the top of the image and east is to the left. 
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Table 1 

Chandra observational data 


ObsID 

(1) 

Date 

(2) 

Instrument 

(3) 

Exposure 

(ks) 

(4) 

Pile-up fraction 

(5) 

Notes 

(6) 


361 

1999-09-20 

ACIS-I 

33.7 

<1% 

on axis 


1302 

1999-09-20 

ACIS-I 

15.7 

<1% 

on axis 


1411* 

1999-10-28 

HRC 

54.0 

- 



378* 

1999-12-30 

ACIS-I 

4.2 

- 

off axis. 

sub-array, X-2 PSF blended 

379* 

2000-03-11 

ACIS-I 

9.1 

- 

off axis. 

sub-array, X-2 PSF blended 

380* 

2000-05-07 

ACIS-I 

5.1 

- 

off axis. 

X-2 PSF blended 

2933 

2002-06-18 

ACIS-S 

18.3 

>10% 

on axis 


6097* 

2005-02-04 

ACIS-S 

58.2 

- 

off axis. 

sub-array, X-2 PSF blended 

5644 

2005-08-17 

ACIS-S 

75.1 

<5% 

on axis. 

sub-array 

6361 

2005-08-18 

ACIS-S 

19.2 

<5% 

on axis. 

sub-array 

8189* 

2007-01-09 

HRC 

61.6 

- 



8505* 

2007-01-12 

HRC 

83.6 

- 



8190* 

2007-06-02 

ACIS-S 

5.8 

- 

off axis. 

X-2 PSF blended 

10027* 

2008-10-04 

ACIS-S 

20.2 

- 

off axis. 

X-2 PSF blended 

10025 

2009-04-17 

ACIS-S 

19.2 

<1% 

off axis. 

sub-array 

10026 

2009-04-29 

ACIS-S 

18.7 

<1% 

off axis 


10542 

2009-06-24 

ACIS-S 

120 

<1% 

on axis 


10543 

2009-07-01 

ACIS-S 

120 

<1% 

on axis 


10925 

2009-07-07 

ACIS-S 

45.1 

<1% 

off axis 


10544 

2009-07-07 

ACIS-S 

74.5 

<1% 

off axis 


11104 

2010-06-17 

ACIS-S 

10.1 

>10% 

on axis 


11800 

2010-07-20 

ACIS-S 

17.1 

<1% 

off axis 


10545 

2010-07-28 

ACIS-S 

96.3 

-5% 

off axis 


13796 

2012-08-09 

ACIS-S 

20.1 

>10% 

on axis 


15616 

2013-02-24 

ACIS-S 

2.1 

<5% 

on axis. 

short observation (2 ks) 

16580 

2014-02-03 

ACIS-S 

47.5 

>10% 

on axis 


17578 

2015-01-16 

ACIS-S 

10.1 

-1% 

off axis. 

sub-array 

16023 

2015-01-20 

ACIS-S 

10.1 

<10% 

on axis 



Note. — Details of the 28 Chandra observations of M82 taken from 1999 to 2015, ordered by date. Column (1) gives the obsID (* 
indicates that this observation was not used in our investigation due to blended PSFs or HRC data), column (2) gives the date of the 
observation, column (3) gives the instrument used, column (4) gives the total exposure time in ks, column (5) gives an estimate of the 
pile-up fraction described in Section!^ and column (6) gives details about if the observation was taken off-axis and/or with a sub-array of 
pixels in order to mitgate the effect m pile-up. We also note if the PSF of X-2 is blended with nearby sources. 
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those of nearby sources. Since our aim is to carry out 
spectral fitting, we do not use the 3 obsIDs of HRC data 
due to the limited spectral capabilities of this instru¬ 
ment. After excluding 6 obsIDs where the PSF of X-2 is 
blended and the 3 obsIDs of HRC data, our dataset for 
spectral extraction consists of 19 ACIS observations. We 
use the CIAO tool specextract to extract the spectra, 
which produces source and background spectra and re¬ 
distribution matrices and auxilliary response files (RMF 
and ARFs). For on-axis observations, we use circular re¬ 
gions of radius 1.2" centered on the source. A radius of 
1.2" encircles 85% of the energy at 6 keV for an on-axis 
point source. This radius increases by a factor of ^ 2 
for off-axis angles of ~ 4', and the PSF elongates into 
an elliptical shape. Therefore for off-axis observations, 
we use elliptical regions with 2" and 1" semi-major and 
semi-minor axes respectively. For observations where X-2 
is faint or absent in the image, we center these regions at 
the same positions relative to the bright persistent source 
X-1 as for the obsIDs where X-2 is clearly detected. We 
also use smaller, 0.8" regions for these observations to 
avoid contamination from the two nearby sources to the 
south. 

In order to carry out background extraction, we ini¬ 
tially considered a large circular region outside the galaxy 
on the same chip as the galaxy, however we find that the 
spectra resulting from the subtraction of this background 



sion. Thus we use the spectrum of the diffuse emission 
for background subtraction. Figure shows the Chan¬ 
dra image from a 75-ks on-axis exposure of M82 (obsID 
5644), where the scaling enhances the diffuse emission. 
This reveals the various structures present in the diffuse 
emission around X-2, which is particularly strong to the 
south east. In our choice of background region, we aim 
to avoid the brightest diffuse emission which is too far 
from the source to contribute to its background, plus to 
avoid the point sources to the south and the PSF of X-1 
to the west. This results in a relatively small 7.6" x 1.9" 
rectangular extraction region to the north east of X-2, 
which though small contains enough signal 20 counts 
ks“^) to undertake the background extraction. We use 
this background region for all obsIDs. 

A major effect that must be taken into account before 
we can perform spectral fitting is pile-up. Pile-up occurs 
when more than one photon lands on the same pixel be¬ 
tween two readouts of the CCD. In this case multiple 
photons are counted as a single event with increased en¬ 
ergy, or they are not counted at all (in the case of grade 
migration). Pile-up therefore affects both the inferred 
flux and spectral shape. The standard readout time for 
the ACIS detectors is 3.2 s, and with a count rate that 
can get up to count s“^, the effect of pile-up on X-2 
observations cannot be neglected. As mentioned above, 
the effect of pile-up can be reduced by the use of a sub¬ 
array of pixels that reduces the readout time to 0.4 s for 
1/8 ACIS-S chip, or by moving the source off-axis, which 
spreads the photons over several pixels. Off-axis obser¬ 
vations have the drawback that the PSF of close sources 
can become blended. 



Figure 2. Figure of obsID 5644 showing the source (circular) and 
background (rectangular) spectral extraction regions. The scaling 
has been modified in this figure with respect to Fig. [^to emphasize 
the diffuse emission. The background region was chosen to contain 
the diffuse emission which was as local to the source as possible 
while also excluding nearby point sources and the PSF of X-1. 
North is up and east is left. 


To estimate the level of pile-up in our data, we use 
the CIAO tool PILEUP_MAP, which outputs an image of 
counts per ACIS frame, which can then be used to esti¬ 
mate a pile-up fractiorQ For off-axis and sub-array ob¬ 
servations, the count rate per frame is typically less than 
0.1, which corresponds to a pile-up fraction of < 5%. 
For on-axis observations with full CCDs, the counts per 
frame are many times higher, with typical pile-up frac¬ 
tions > 10%. We list these estimates in Table (T] We 
identify 4 observations where the estimated pile-up frac¬ 
tion is > 10%, 2933, 11104, 13796 and 16580. The ef¬ 
fect of pile-up can be accounted for in spectral fitting for 
cases where pile-up is not too strong, which we judge to 
be the case here, using a convolutio n model based on an 
algorithm presented in Davis (2001). 

We use XSPEC (vl2.k2) to carry out spectral fitting 
of X-2. We group the spectra with a minumum of 20 
counts per bin using the HEASARC tool GRPPHA using 
the fif statistic and carrying out background subtrac¬ 
tion. There are not enough counts in the spectra ex¬ 
tracted from obsID 1302 and 15616 for fitting, thus 
we group the spectrum to have a minimum of 1 count per 
bin and use the Cash fit statistic. We fit the spectra in 
the energy range 0.5—8 keV with two models, a power- 
law (powerlaw) model and a disk black body (diskbb) 
model, both of which are subjected to photoelectric ab¬ 
sorption local to the source (zwabs). 

Despite subtracting off the diffuse background close to 
the position of X-2, we find that soft excess emission 
still persists over the above models. It is possible that 
this is residual diffuse emission not accounted for by the 
background subtraction, however it may also originate 


^ http://cxc.harvard.edu/csc/memos/files/Davis_pileup.pdf 
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from optically-thin plasma that has been photo-ionized 
by the pulsar, evidence for which has been found around 
other ULXs such as Holmberg II X-1 (?), NGC 7424 
ULX 2 (?) and NGC 5408 X-1 (?). We account for 
this soft excess component in any case by adding the 
collisionally-ionized diffuse gas model, apec. While many 
X-ray pulsars indeed present soft X-ray excesses thought 
to be produced by the reprocessing of the hard X -rays in 
the inner r egion of the accretion disk (see e.g. Hickox 
1004), the absorption in this system (A^h^ 3 x 10^^ 
makes it unlikely that the excess we observe is 
intrinsic to X-2. Similarly, we find that if we subject 
the apec model to the same absorption that is applied 
to the power law or diskbb model, it cannot account for 
the soft excess. Furthermore, adding absorption to this 
component that is independent from the absorption of 
the pulsar does not produce an improvement in We 
conclude that the soft excess most likely originates from 
foreground material in the line-of-sight that is sufficiently 
outside the bulk of the diffuse gas that it is not absorbed 

The absorption column and F (or Tin) are degener¬ 
ate given a limited bandpass (i.e. a hard spectrum can 
be fitted with a steep F and high A^h or a fiat F and 
low A^h)- For this reason we simultaneously fit for A^h 
across all epochs, which is thus driven by the spectra with 
the greatest counts. The parameters of the powerlaw or 
diskbb model are free for each epoch. The absorption 
in the spectrum is attributed to the interstellar medium 
in M82 rather than being local to the neutron star, and 
thus is not expected to change on the time-scales that we 
are considering. We confirm this explicitly by fitting the 
spectra individually with free A^h parameters, and find 
no evidence for a variable absorber. Furthermore, as we 
attribute the soft excess to diffuse emission, we also si¬ 
multaneously fit for the temperature of the apec model. 
We leave the normalizations of this model free since the 
small circular source extraction regions for on-axis ob¬ 
servations and the larger elliptical extraction regions for 
off-axis observations contain differing amounts of the dif¬ 
fuse emission. 

For both model cases we test the effect of adding the 
multiplicative pileup model. In these cases, the frame 
time parameter is set dependent on the size of the sub¬ 
array used. We set the PSF fraction (not the fraction of 
the PSF included in the extraction region but the fraction 
of counts in the region which are from the point source 
whose pile-up is being modeled) to 95%. The only pa¬ 
rameter left free when using this model is a, which is the 
grade morphing parameter. The parameter a is related 
to the grade migration function, G = where p is the 

number of piled photons. We consider four factors when 
deciding whether to include this model in the spectral fit. 
We visually inspect the spectrum for a telltale turn-up at 
high energies. We use the estimated pile-up fraction from 
the PILEUP_MAP tool and consider if the observation was 
taken off-axis and/or with a sub-array. Generally on-axis 
observations show pile-up fractions greater than 10% and 
we thus include the pileup model in the fit. Lastly, if the 
pileup model is included and the best-fit a parameter 
converges on a small number or zero, which is considered 
unphysical, we remove the pileup model from the fit. 

For the powerlaw model, the free parameters are the 
photon-index, F, and the normalization. We allow F to 


vary between -3 and 10. For the diskbb model, the free 
parameters are Tin, the temperature of the inner disk in 
keV, and the normalization. We restrict the disk temper¬ 
ature to < 10 keV. For both models the normalization of 
the apec model is a free parameter. As described above, 
the temperature of this component we fit for simultane¬ 
ously across all epochs and found to be keV 

and 0.58^0 17 heV for the powerlaw and diskbb models 
respectively. The A^h of the zwabs model is also fitted 
simultaneously and found to be x 10^^ cm“^ and 

2.8lo‘j9 ^ cm“^ respectively. The diskbb model re¬ 
quires less absorption since it predicts fewer soft X-ray 
photons than the powerlaw model. The redshift is fixed 
at 0.00067 for all model components. 

Uncertainties on the free parameters are calculated us¬ 
ing the Ax^=4.61 criterion, which corresponds to 90% 
confidence level for two interesting parameters. 


3. NUSTAR DATA AND ANALYSIS 


NuSTAR ( Harrison et al.|20r3 ) observed the M82 field 
7 times between 2014 January 23 and 2014 March 06, as 
described in B14, for a total exposure of 1.91 Ms during 
which the pulsations were detected. A Chandra expo¬ 
sure (47.5 ks, obsID 16580) overlapped with the NuS¬ 
TAR observation during a period where the pulsations 
were present. While B14 presented timing and photo¬ 
metric analysis of the pulsar, we aim here to present for 
the first time some of its spectral characteristics above 10 
keV by isolating the pulsed component. For this we used 
the NuSTAR data analysis software (NuSTARDAS) ver¬ 
sion 1.2.0 and NuSTAR CALDB version 20130509 with 
the standard filters to obtain good time intervals, ex¬ 
cluding the periods where the source was occulted by 
the Earth or was transiting through the South Atlantic 
Anomaly. We used the pulsar ephemeris described in 
B14 to extract “pulse-on” and “pulse-off” spectra. The 
pulse-on spectrum is defined to be the brightest 25% of 
the pule profile, while the pulse-off is defined to be the 
faintest 25% of the pulse profile. We then subtract the 
pulse-off spectrum from the pulse-on spectrum to obtain 
the pulsed spectrum which we model with some simple 
models in order to characterize the data and to facilitate 
comparison with other well-studied pulsars. Due to the 
triggered read out of the NuSTAR detectors, pile-up is 
not an issue at the flux levels of M82. The NuSTAR 
data were rebinned to a signal-to-noise of 3, providing a 
significant signal up to ^ 40—50 keV. We fit the pulsed 
spectrum with models consisting of an absorbed power- 
law continuum both with and without an exponential 
cut-off (cutoff pi), both of which are subjected to inter¬ 
stellar absorption, where we fix the Ar to 3 x 10^^ cm“^, 
as determined from the Chandra analysis. 


4. RESULTS 

The results of the spectral fits using both the power- 
law and disk black-body model are given in Table The 
fit statistic, for the combined fits with the power-law 
model is 1439.89 with 1315 degrees of freedom (x^=1.09), 
whereas for the disk black body model =1323.03 with 
1315 degrees of freedom (y^=1.01). Comparing these 
fit statistics suggests that the disk black body model is 
overall the best model. However, at the lowest luminosi¬ 
ties (^ 10^^ ergs“^), the temperature of the disk black 
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body model is not well constrained and hits the upper 
bound of 10 keV imposed for spectral fitting during the 
error calculations. In some cases the error calculations 
fail. This is most likely due to the low number of counts 
at these luminosities. 

The intrinsic luminosity estimates (corrected for ab¬ 
sorption) between the two models differ due to both the 
diverging spectral shape above 8 keV (luminosity mea¬ 
surements are extrapolated to 10 keV for comparison 
with previous works) and in the soft X-ray band, plus 
the differing pile-up estimates caused by the difference 
in spectral shape of the models. The disk black body 
model gives a systematically lower 0.5—10 keV intrinsic 
luminosity, by typically 25%. Nonetheless, our analysis 
shows that X-2 is frequently observed to be emitting well 
above its Eddington luminosity assuming a typical NS 
mass of 1.4 Mq, regardless of the model used. This is il¬ 
lustrated in Figure which shows intrinsic (unabsorbed) 
Lx against time, for 9/19 (47%) observations that we 
analyzed, we found that X-2 emits at a luminosity in ex¬ 
cess of 10 ^^ ergs“^ and thus relatively persistent, ra ther 
than transient as i dentified by Feng & Kaaret (2007) and 
Kong et al. (2007). 


We also note the period of extreme flux variability in 
2010 where the source drops almost two orders of magni¬ 
tude in brightness in the space of a month. The exposure 
times of these two observations (11104 and 11800) are 
too short to reveal any significant variability during the 
observations, however. Of the longer observations avail¬ 
able, obsIDs 5644 (75 ks) and 10545 (96 ks) are least 
affected by pile-up due to a sub array of pixels used in 
the former and X-2 lying at off-axis angles in the latter. 
The light curves of these observations are shown in Fig¬ 
ure]^ While X -2 does not show much variability during 
obsID 5644, it does show variability of up to a factor of 
two on ks-timescales during obsID 10545. 

We also examine the spectra of X -2 for these two long 
observations in order to gain further insights into the 
source, shown in Figure and Figure During obsID 
5644 X-2 was observed to be near its peak luminosity 
of ^ 10^^ ergs“^, while during obsID 10545 X-2 was at 
a lower luminosity of ^ 4 x 10^^ ergs“^. In Figure 
we show the data-to-model ratios of these spectra when 
fitted by the power-law and disk black body models. In 
neither case do the ratios show any deviations indicative 
of a bad fit. Furthermore, the values are comparable 
between the models. 

Figure shows the relationship between the spectral 
parameters T of the power-law model and Tin of the 
disk black body model with respect Tx- For Lx> 10 ^^ 
ergs“^ the mean T = 1.33 ±0.15 (Icr), whereas the mean 
Tin = 3.24 ±0.65. Below 10^^ ergs“^, T and Tin are not 
well constrained and show a large spread in values. 

In addition to the Chandra spectral analysis, we have 
conducted NuSTAR pulse-phased spectroscopy of the 
pulsar in the 3—50 keV range, as described in Section ^ 
We fit the pulsed spectrum with a power-law model anda 
model with an exponential cut-off (cut off pi in XSPEC), 
both of which are subjected to photo-electric absorption. 
We find that the power-law with a cut-off is significantly 
preferred (Ax^= 20 for one additional free parameter). 
The fit is excellent (y^/DoF = 132/126) and we find the 
photon index of the pulsed component to be T = 0.6±0.3, 


with a high energy cut-off, Eq = 14^3 keV. The average 
pulsed flux in this band is 5.7±0.4 x 10“^^ erg cm“^ s“^, 
corresponding to a luminosity of 7.5 x 10^^ ergs“^ at 3.3 
Mpc. We present the pulsed spectrum of M82 X -2 in 
Figi unfolded through the instrumental response as¬ 
suming the cut-off power-law model. 

5. COMPARISON OF THE M82 PULSAR WITH OTHER 
LUMINOUS X-RAY PULSARS AND IMPLICATIONS FOR 
THEORETICAL MODELS 

Other examples of luminous (> 10^^ ergs“^) X-ray 
pulsars include SMC X- 1 , LMC X-4, GRO J1744-28, 
RX J0059.2-7138, XTE JOlll.2-7317 and A0538-66. Al¬ 
though M82 X-2 can reach luminosities an order of mag¬ 
nitude brighter than these sources, a comparison with 
them is valuable. The X-ray spectra of these sources are 
typicaly fitted with a power-law model with T = 0.5 —1.5, 
subjected to absorption along the line of sight, and in 
some cases with a high-energy cut-off with energies rang- 
i ng from 5—30 keV. 

Paul et al.| (| 2002 |) fit the phase-averaged spectra of 
[ LIVIC X-4 in the 0.7—10 keV band using 
The spectrum of SMC X -1 was modeled 


SMC X-1 and 
ASCA data 

well with a cut-off power-law, where T = 0.91 ± 0.03 and 
the cut-off energy is 5 . 5^0 5 keV. The spectrum of LMC 
X-4 did not require a cut-off, and could be reproduced 
by a plain power-law with T = 0.69 ±0.04. The other ex¬ 
amples of super-Eddington pulsars, GRO J1744-28, RX 
J0059.2-7138, XTE JOlll.2-7317 and A0538-66 are also 
well described by cut-off power-laws or power-law mod¬ 
els, both with similar parameters (see Nishiuchi et al. 
19991 lYounes et al.||2015[ ISidoli et al.||2015[ |Yokoga ~ 


et al. ]|2000[ |Skinner et alj|1982[ repectively). 


T'he phase-averaged 0.5—8 keV X-ray spectral proper¬ 
ties of M82 X-2 are very similar to these other sources, 
with r = 1.33 ±0.15 for a power-law model with no cut¬ 
off. The long Chandra observation obsID 5644 where the 
source is caught in an ultra-luminous state, but where 
pile-up is negligible, offers us the opportunity to test 
for the presence of a cut-off and if any constraints can 
be placed on it. Fitting this spectrum with this model 
yields T = 0.70lo‘65 = 6 .I 9 I 2 V keV, where 

=389.87 for 338 degrees of freedom. The fit with the 
power-law without a cut-off gives x^ =395.66 for 339 de¬ 
grees of freedom, thus the addition of the cut-off yields 
an improvement in x^ of only 6 for the addition of one 
parameter and the cut-off energy is not well constrained 
in the Chandra data alone. 

NuSTAR^ however, can measure the cut-off due to its 
sensitivity above 10 keV. From analysis of the pulsed 
component in the 3—50 keV range we find that the pulsed 
spectrum is best fitted by a power-law with a high-energy 
cut-off, where T = 0.6 ± 0.3 and Eq = 14^3 keV. These 
values are very similar to the phase-averaged spectrum 
described above, albeit with much better constraints ow¬ 
ing to the high-energy sensitivity of NuSTAR. The sim¬ 
ilarity with the phase-averaged spectrum from Chandra 
indicates that the pulsed spectrum dominates the emis¬ 
sion of the pulsar. It is also interesting to note that the 
cut-off energy of the pulsed component is higher than 
that observed in other ULXs where the nature of the ac- 
cretors remain unknown, whose spectra typically cut off 

at 6-8 keV (e.g. iBachetti et al.||2013l iWalton et al.||2013 

. " - 
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Table 2 

Spectral fitting results 


ObsID 

Date 

Power-law model parameters 

Disk black body model parameters 



Qipl 

Tpl 

Ux,pi 

<Odisk 

^in,disk 

Ux,disk 





(xlO^® ergs“^) 


(keV) 

(xlO^® ergs“^) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

361 

1999-09-20 

- 

TsoTTir 

1 n+ 0.6 

-L- 0 _q 4 

- 

2-64115S 

o.sl^i 

1302 

1999-09-20 

- 

1 '71 + 1.53 

l-+0°5 

- 

- 

- 

2933 

2002-06-18 

0.691^^9 

1 [;:p:+0.21 

i.OO_o 20 



2-46lg;Il 

579lif 

5644 

2005-08-17 

- 

1 42+0-05 
-L-^^_0.05 

1151:1 

- 

2.86lol° 

921^ 

6361 

2005-08-18 

- 

1 04 +O.II 

10211 

- 

3.16l°l® 

83l® 

10025 

2009-04-17 

- 

2 41+0-90 

7 q + 7.3 
'•'^-3.0 

- 

1-781S.84 

3 6 “^^'® 
'^•'^- 1.0 

10026 

2009-04-29 

- 


r, 0 + 13.6 
^•'^- 1.6 

- 

- 

- 

10542 

2009-06-24 

- 

1 QQ+0-28 

i .»»_0 28 

q c:+0.5 
0-0-04 

- 

1 94+0-67 
-*-•^^-0.40 

9 q+ 0.1 

^•'^- 0.2 

10543 

2009-07-01 

- 

i.ZU_Q g 4 

1 4+0-2 
-'■•^- 0.2 

- 

5 91 +^ 
^•^^-3.03 

1 9 + 0.3 
^•^-0.3 

10925 

2009-07-07 

- 

1 08+0-08 
i.08_o 76 

^•°-0.5 

- 

7.42l“ 

9 q+0.5 
^•'^-0.7 

10544 

2009-07-07 

- 


1 1 +°-^ 
^•^-0.3 

- 

9 .O 1 II 

1 1+0.3 
^•-^-0.4 

11104 

2010-06-17 

0.2o1°;29 

1 14+0.29 

197+29 

-64 

n 97 + 0.33 
-0.24 

4-23ll4g 

143+180 

1^0-40 

11800 

2010-07-20 

- 

0 c:-] +1-25 

4.112.1 

- 

l-SSlS.gg 

-^•^- 0.8 

10545 

2010-07-28 

- 

1 34+0-09 
-'-•'^^-0.09 

36li 

- 

q gg + 0.51 

o.^^-O 38 

soli 

13796 

2012-08-09 


i-i+r7o 

227l°5 

n 25“*”®'^^ 

0 '-^O_Q Qg 

4 17+1-38 
-0.83 

1971®® 

15616 

2013-02-24 

- 



- 

- 

- 

16580 

2014-02-03 

-0.06 

-| Q7+0.09 
-0.09 

997+61 

-49 

0'-^' -0.05 

q qq + 0.70 
O.OO-o 51 

2101^7 

17578 

2015-01-16 

- 

1 ' 55 + 0-22 
i.OO-o 23 

5211 

- 

2 3«+0-78 
^•OD -0 45 

39li 

16023 

2015-01-20 

- 

1 90 + 0.27 
-^•^^-0.28 

76i^r 

- 

3.3711®® 

36l“ 


Note. — Results of the X-ray spectral fitting of the 19 obsIDs where reliable spectral information on X-2 could be extracted, fitted with 
an absorbed power-law and absorbed disk black body models. Column (1) gives the obsID, column (2) gives the date of the observation, 
column (3) gives the grade morphing parameter of the pileup model when convolved with the power-law model, if this is used in the fit 
(‘-’ indicates the pileup model was not used), column (4) gives the photon index of the power-law model, column (5) gives the unabsorbed 
0.5—10 keV luminosity of the power-law model. Column (6) gives the grade morphing parameter of the pileup model, when convolved with 
the disk black body model, if this is used in the fit (‘-’ indicates the pileup model was not used), column (7) gives the disk temperature in 
keV. ‘n’ indicates that this parameter hit the upper limit of 10 in the spectral fit and ‘T indicates it hit the lower limit of 0. indicates 
where the error calculation fails. Column (8) gives the unabsorbed 0.5—10 keV luminosity of the disk black body model. 



papers that have aimed to explain the nature of the ultra- 
luminous pulsar have presented different scenarios for the 
strength of its magnetic field. Since E^yc is directly pro¬ 
portional to the strength of the magnetic field, then Eq 
could potentially yield information about the magnetic 
field strength. Following this, the 14 keV cut-off that we 
measure in the pulsed spectrum implies a 5 ^ 10^^ G 
magnetic field. However, we note that there are excep¬ 
tions to the above relation, for exam ple KS 1947+319 , 
wher e £'cvc=12 keV and F^c=22 keV ( Furst et al. 12014 ). 


DairOsso et al. (|2015) also discuss the variability ex¬ 
hibited by X-2, which can be explained by relatively 
small changes in the mass accretion rate in the pres¬ 
ence of a strong magnetic field, whereby the source at 
low luminosities enters the propeller regime. We note, 
however, that the minimum luminosity at which accre¬ 
tion is possible in the presence of a H = 10^^ G magnetic 
field, a strength which t hey favor, with a 1 .37 s period is 
^ 2 X 10^^ ergs“^ (see Stella et ^|1986). Gonsidering 


that we observe X-2 at much lower luminosities, argues 
against such a strong magnetic field given the discussion 
above. 

Goncerning the duty cycle of X-2, we have found that 


the source is more per siste nt than prev i ously reported by 
Feng & Kaaret|(2007) and Kong et al. (2007), due to the 
longer baseline of our investigation and larger data set. 
For 9/19 (47%) observations that we analyzed, we found 
that X-2 emits at a luminosity in excess of 10^^ ergs“^. 
Luminosities of ^ 10^^ — 10^^ ergs“^ imply that the 
neutron star is growing at a rate of ^ 2 x 10“^ — 10“^ 
Mq yr“^, assuming isotropic emission and a mass-to- 
energy conversion efficiency of unity, meaning it will col¬ 
lapse into a black hole within ^10 — 100 million years. 
These results could have important implications for the 
formation and growth of supermassive black holes, theo¬ 
retical modeling of which often employs an early super¬ 
critical accretion phase to explain the masses of the 
super-massive black holes found in quasars a t z ^ 1 (e.g. 
Volonteri fc Rees||2QQ6 Mortlock et al.||2Qll ). Low-mass 
X-ray binaries have also been proposed to be a potential 
source of ionizing radiation for heating the intergalactic 
mediu m during the epoch of reionization (Fragos et al. 
2013), which may indeed have a contribution from a non- 
negligible population of ultra-luminous pulsars. 

6. SUMMARY AND CONCLUSIONS 


In this paper we have conducted a temporal and spec¬ 
tral analysis of the 0.5-8 keV X-ray emission from the 
ultra-luminous X-ray pulsar M82 X-2 from 15 years of 
Chandra observations and pulse-phased spectroscopy in 
the 3—50 keV band from NuSTAR data. Our main find¬ 
ings are as follows: 
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Figure 3. Long-term activity of X-2 over the 15 years of Chandra observations showing that the source is frequently observed to be 
radiating at many times its Eddington limit. Vertical lines indicate the 90% confidence range on the measured 0.5—10 keV luminosity, 
which is calculated from the power-law model assuming a distance to M82 of 3.3 Mpc. Dark blue squares show observations which were 
taken off-axis and/or with a sub-array of pixels to mitigate the effect of pile-up. Light blue squares are observations taken on-axis with the 
full array of CCDs. The horizontal lines show the Eddington limiting luminosity for a 1.4 Mq object, along with 10 and 100 times this. 


0.20 

cn 

^ 0.15 

0) 

cd 

i 0.10 

J.,., 


0.00 

0 20 40 60 80 100 

Time (ks) 




5644, 2010-07-28 


10545, 2005-08-17 

3 


11 igi 




Figure 4. Light-curves of X-2 for obslDs 5644 and 10545, with 2 
ks bins. 


• When fitted with a power-law model, the average 
photon index for epochs where Lx> 10^^ ergs“^ 
in the 0.5—8 keV band is F = 1.33 ± 0.15. For the 
disk black body model, the average temperature is 
Tin = 3.24 ± 0.65 keV. This spectral shape is con¬ 
sistent with other luminous X-ray pulsars. We also 
investigated the inclusion of a soft excess compo¬ 
nent and spectral break finding that the spectra 
are also consistent with these features common to 
luminous X-ray pulsars. 

• The pulsed emission of X-2 in the 3—50 keV band 
from NuSTAR data is best fitted by a power-law 
with a high-energy cut-off, where F = 0.6 ±0.3 and 
Ec = 14^3 keV with a luminosity of 7.5 x 10^^ 
ergs“^. 
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Figure 5. 75 ks Chandra spectrum of X-2 taken in 2005 (obsID 
5644), taken on axis but with a sub-array of pixels to reduce pileup. 
The top panel shows the unfolded spectrum fitted with an absorbed 
power-law, plotted with a dashed line. The dotted line shows the 
apec model used to model the excess diffuse background. The 
middle panel shows the data-to-model ratio of this fit. The bottom 
panel shows the data-to-model ratio of a fit with the disk black 
body model. 
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Figure 6. 96 ks Chandra spectrum of X-2 taken in 2010 (obsID 
10545), taken off axis to reduce pileup. The top panel shows the 
unfolded spectrum fitted with an absorbed power-law, plotted with 
a dashed line. The dotted line shows the apec model used to model 
the excess diffuse background. The middle panel shows the data-to- 
model ratio of this fit. The bottom panel shows the data-to-model 
ratio of a fit with the disk black body model. 

• Our results show that X-2 has been remarkably 
active over the 15-year period considered. We 
find that for 9/19 (47%) observations that we an¬ 
alyzed, the pulsar appears to be emitting at a lu¬ 
minosity in excess of 10^^ ergs“^, which is greater 
than 10 times its Eddington limit. Luminosities of 
^ 10^^ — 10^^ ergs“^ imply that the neutron star 
is growing at a rate of ^ 2 x 10“^ — 10“^ Mq yr“^ 
and is expected to collapse into a black hole within 
^10 — 100 million years. 
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Figure 7. Figures showing (top) the relationship between E and 
Lx of X-2 when fitted with a power-law model and (bottom) the 
relationship between Tin and Lx when fitted with a disk black body 
model. Dark blue squares show observations which were taken off- 
axis and/or with a sub-array of pixels to mitigate the effect of 
pile-up. Light blue squares are observations taken on-axis with the 
full array of CCDs. In some cases, at the lowest luminosities, the 
best-fit disk temperature is off the scale of this figure and with 
large unconstrained values. These are indicated with vertical lines 
in the figure not associated with a square. 


2013)/ERC Grant Agreement n. 617001. 
Facilities: Chandra (ACIS), NuSTAR 

REFERENCES 


Bachetti, M., Rana, V., Walton, D. J., et al. 2013, ApJ, 778, 163 
Bachetti, M., Harrison, F. A., Walton, D. J., et al. 2014, Nature, 
514, 202 

Basko, M. M., &: Sunyaev, R. A. 1975, A&A, 42, 311 
—. 1976, MNRAS, 175, 395 

Christodoulou, D. M., Laycock, S. G. T., & Kazanas, D. 2014, 
ArXiv e-prints, arXiv:1411.5434 
DairOsso, S., Perna, R., & Stella, L. 2015, MNRAS, 449, 2144 
Davis, J. E. 2001, ApJ, 562, 575 

Ek§i, K. Y., Andag, i. C., Qikmtoglu, S., et al. 2015, MNRAS, 
448, L40 

Endo, T., Nagase, F., & Mihara, T. 2000, PASJ, 52, 223 
Feng, H., & Kaaret, P. 2007, ApJ, 668, 941 
Foley, R. J., Fox, O. D., McCully, C., et al. 2014, MNRAS, 443, 
2887 

Fragos, T., Lehmer, B. D., Naoz, S., Zezas, A., & Basu-Zych, A. 
2013, ApJ, 776, L31 

Fragos, T., Linden, T., Kalogera, V., & Sklias, P. 2015, ApJ, 802, 
L5 

Fiirst, F., Grefenstette, B. W., Staubert, R., et al. 2013, ApJ, 

779, 69 

Fiirst, F., Pottschmidt, K., Wilms, J., et al. 2014, ApJ, 784, L40 
Gandhi, R, Isobe, N., Birkinshaw, M., et al. 2011, PASJ, 63, 505 




























































Chandra and NuSTAR observations of the ULX pulsar in M82. 


11 



Energy (keV) 


Figure 8. The NuSTAR 3—50 keV unfolded pulsed spectrum of 
M82 X-2, fitted with a power-law, E = 0.6±0.3, with a high energy 
cut-off, Eq = 14^3 keV. The average pulsed flux in this band is 

5.7 ± 0.4 X 10“^^ ergcm“^ s“^, corresponding to a luminosity of 
7.5 X 10^^ ergs“^ at 3.3 Mpc. 


Harrison, F. A., Craig, W. W., Christensen, F. E., Hailey, C. J., 
& Zhang, W. W. 2013, ApJ, 770, 103 
Hickox, R. C., Narayan, R., & Kallman, T. R. 2004, ApJ, 614, 881 
Kaaret, P., Simet, M. G., & Lang, C. C. 2006, ApJ, 646, 174 


Kluzniak, W., & Lasota, J.-P. 2015, MNRAS, 448, L43 
Kong, A. K. H., Yang, Y. J., Hsieh, P.-Y., Mak, D. S. Y., & Pun, 
C. S. J. 2007, ApJ, 671, 349 
Lyutikov, M. 2014, ArXiv e-prints, arXiv: 1410.8745 
Makishima, K., Mihara, T., Ishida, M., et al. 1990, ApJ, 365, L59 
Matsumoto, H., Tsuru, T. G., Koyama, K., et al. 2001, ApJ, 547, 
L25 

McDonald, A. R., Muxlow, T. W. B., Wills, K. A., Pedlar, A., & 
Beswick, R. J. 2002, MNRAS, 334, 912 
Mortlock, D. J., Warren, S. J., Venemans, B. P., et al. 2011, 
Nature, 474, 616 

Mukherjee, E. S., Walton, D. J., Bachetti, M., et al. 2015, ArXiv 
e-prints, arXiv: 1502.01764 

Mushtukov, A. A., Suleimanov, V. F., Tsygankov, S. S., & 
Poutanen, J. 2015, ArXiv e-prints, arXiv: 1506.03600 
Nishiuchi, M., Koyama, K., Maeda, Y., et al. 1999, ApJ, 517, 436 
Pasham, D. R., Strohmayer, T. E., & Mushotzky, R. F. 2014, 
Nature, 513, 74 

Paul, B., Nagase, F., Endo, T., et al. 2002, ApJ, 579, 411 
Rana, V., Harrison, F. A., Bachetti, M., et al. 2015, ApJ, 799, 121 
Ranalli, R, Comastri, A., Origlia, L., & Maiolino, R. 2008, 
MNRAS, 386, 1464 

Shao, Y., & Li, X.-D. 2015, ApJ, 802, 131 
Sidoli, L., Palombara, N. L., Esposito, P., Tiengo, A., & 
Mereghetti, S. 2015, MNRAS, 449, 3710 
Skinner, G. K., Bedford, D. K., Eisner, R. F., et al. 1982, Nature, 
297, 568 

Stella, L., White, N. E., & Rosner, R. 1986, ApJ, 308, 669 
Volonteri, M., & Rees, M. J. 2006, ApJ, 650, 669 
Walton, D. J., Fuerst, F., Harrison, F., et al. 2013, ApJ, 779, 148 
Weliachew, L., Fomalont, E. B., & Greisen, E. W. 1984, A&A, 
137, 335 

Yokogawa, J., Paul, B., Ozaki, M., et al. 2000, ApJ, 539, 191 
Younes, G., Kouveliotou, C., Grefenstette, B. W., et al. 2015, 
ApJ, 804, 43 





























